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Abstract. We present a new algorithm to search for dis- 
tant clusters of galaxies on catalogues deriving from imag- 
ing data, as those of the ESO Imaging Survey. 

Our algorithm is a matched filter one, similar to that 
adopted by Postman et al. 1996, aiming at identifying clus- 
ter candidates by using positional and photometric data 
simultaneously. The main novelty of our approach is that 
spatial and luminosity filter are run separately on the cat- 
alogue and no assumption is made on the typical size nor 
on the typical M* for clusters, as these parameters inter- 
vene in our algorithm as typical angular scale a and typical 
apparent magnitude m* . Moreover we estimate the back- 
ground locally for each candidate, allowing us to overcome 
the hazards of inhomogeneous datasets. As a consequence 
our algorithm has a lower contamination rate - without 
loss of completeness - in comparison to other techniques, 
as tested through extensive simulations. We provide cat- 
alogues of galaxy cluster candidates as the result of ap- 
plying our algorithm to the I-band data of the ElS-wide 
patches A and B. 

Key words: Methods: numerical - Galaxies: clusters: 
general - Cosmology: large-scale structure of Universe 



1. Introduction 

The quest for high-redshift [z 0.5) clusters of galaxies 
has recently received a lot of attention, and several search 
methods have been put forward. The mere existence of 
rich clusters of galaxies at very high rcdshifts is cosmo- 
logically relevant and their importance as a discriminant 
among different theoretical cosmological models is widely 
ackn owledged (see eg. Via na fc L iddle 1996| ; Carlber g et a l. 



19971 ; Bahcall, Fan & Cen |1997[ Bartelmann et al. |1998| ) . 



On the other hand, a detailed study of the properties of 
clusters, as well as of their member galaxies, at high red- 
shifts and the confrontation of these with local well-known 
systems gives precious insight on their formation and on 
the evolution of their properties with redshift and envi- 



ronment (Dressier et al. 1997; van Dokkum et al. 1998a, 
1998b[ S mail et al. |1998] ; Morris et al. |1998t Stanford et 
al. 1998 and references therein). The now well-established 



evolution of galaxies with redshift still lacks a clear under- 
standing of the physical processes that guide it: clusters 
are privileged observational targets to distinguish between 
intrinsic and environmental effects. 

To perform the search for distant clusters of galaxies, 
various authors have been taking advantage of several win- 
dows in the spectral emission of the different cluster com- 
ponents, namely galaxies and gas, and applying the best 
suited method in each case. Thus, optical, near IR, X-ray 
and even variations in the cosmic microwave background 
emission have been used to detect cluster candidates at 
distant rcdshifts. All different methods have advantages 
and drawbacks and the key is to view them as comple- 
mentary. Describing each one of them and presenting their 
respective pros and cons is beyond the scope of this paper. 
Nevertheless, it is useful to recall what has been done up 
to now, in order to place our work in the actual scenario 
of these researches. 

Among the different techniques that basically search 
for surface density enhancements on the galaxy 2D spa- 
tial distribution - adding or not magnitude information as 
well -, some of the most popular ones are the application 
of adapted filters (eg. kernel, wavelets, matched filters), 
counts-in-cells techniques, percolation algorithms or even, 
lately, Voronoi tessellation, on optical photometric data 
(eg. Schectman 1985| ; Lumsden et al. 1992; Dalton et al. 
1994|; Escalera & MacGillivray |1994 |199€ ; Pisani 



Lidman & Peterson |1996 
Gal et al. 



al. 1996 



2000 



Fadda et al. 
Ostrander et al. 



NIR data (Stanford et 



1996 



1998; Postman et 



1998), and on 



1998, Ramclla, Nonino , 



1997 ; Mendes de Oliveira et al. 
Boschin 1998[ ) . Apart from these 



more or less elaborate techniques, simpler high-contrast 
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methods still prove to be successful (Couch et al. 1991), 
but they probably don't sample adequately the full dis- 
tant cluster population and their selection criteria are not 
well defined. Some of these algorithms benefit also from 
- or are mainly based upon - multiband colour informa- 
tion that helps i solati ng red galaxies at higher redshifts 
(Gladders & Yee |2000|) . 
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Very significant results have been achieved via the 
search for sources of extended emission in X-ray surveys 
with detection algorithms which are designed to probe a 
broad range of cluster parameters such as X-ray flux, sur- 



face brightness and morphology (eg. Henry et al. 1992 with 
EMSS; RIXOS by Castander et al. [l995] the RDCS by 



Rosati et al. 1995, 1998; NEP by Henry et al. 1997; Gioia 



1997; WARPS by Scharf et al. 1997; Jones, Scharf, Ebel 



ing et al. |1997t SHA RC by Collins et al. |1997| the BCS by 
Ebeling et al. [1998 ; the Cf A large area survey by V ikhlin in 



et al. |1998at |1998b| ; REFLEX by Bohringer et al. |199J). 

Other emerging strategies for cluster search include: 
the detection of extragalactic background light fluc- 
tuations in shallow optical images (Dalcanton 1996; 
Zaritsky et al. 1997); deep imaging around privileged 



sites of density enhancements such as distant powerful 
radio-galaxies or radio-loud quasars (eg. Le Fevre et 
al. 1996; Deltorn et al. 1997); narrow-band imaging 



to search for concentrations of Ly-a emitters around 
previously known weak radio QSO's at high redshift (eg. 
Pascarelle et al. 1996; Campos et al. 1999). Recently, 
the detection of decrements in the cosmic microwave 
background radio emission were also attributed to the 
presence of distant gaseous syste ms, possibly c luster s 
(Jones, Saunders, Baker et al. 1997 ; Richards et al. 1997 ), 
that scatter the microwave background radiation via th e 
Sunyaev-Zel'dovich effect (Sunyaev & Zel'dovich 1980|) . 
One should mention, though, that these strategies for 
cluster search refer to z > 2 systems and we ignore 
"whether these are massive, collapsed systems, groupings 
within unvirialized "sheets" of galaxies, or collections of 
"protogalactic" fragments destined to merge into single, 
more massive galaxies" (Dickinson 1996| ). 



If one aims at producing statistically significant re- 
sults, a good catalogue of clusters, preferably spanning a 
large interval of redshift, obtained using well defined selec- 
tion criteria, over a reasonable sky area, is needed. Our al- 
gorithm, by uniformly detecting clusters over a wide range 
of redshifts and cluster sizes, is well suited to provide such 
a catalogue. We stress that it does not aim, however, at 
performing accurate estimates of neither the redshift nor 
the richness of the cluster candidates. 

The ESO Imaging Survey (EIS, Renzini & da Costa 



1997), covering a final area of 17 square degrees of the 



southern sky in the I-band, up to limiting magnitude 
I ~ 23, provided us with a good opportunity to produce 
such a catalogue. We have developed an automated cluster 
search algorithm and applied it to the catalogues derived 
from the EIS imaging data to obtain a reliable set of clus- 
ter candidates up to estimated z ~ 1.1. Having a list of 
robust candidates is highly desirable, before proceeding to 
the spectroscopic observations with very large telescopes, 
and our algorithm has a high success rate - that is, a high 
completeness level - without being overwhelmed by con- 
tamination in the form of spurious detections. We have 



already selected three of the highest redshift candidates 
(z > 0.5) and performed the spectroscopic follow-up with 
VLT, confirming their real existence as physically bound 
systems (discarding the possibility of false chance align- 
ments) and determined their distance/redshift with accu- 
racy. 

In this paper we describe our algorithm in section || 
stressing its strong points and advantages relatively to 
others present in the literature. Applying it to the EIS- 
wide data of patches A and B produced a catalogue of 
candidate clusters of galaxies that we present in section [s| 
This set is quite different from the one produced by the 
ElS-team of ESO for the same original data set but using 
a different approach and in section |I] we fully investigate 
all possibilities that may account for such a discrepancy. 
Final remarks are given in section |^. 

We shall use H = 50 km s _1 Mpc _1 and go = 0.5 
throughout this paper, unless explicitly stated otherwise. 

2. Our new algorithm 

In their pioneering work on the Palomar Distant Clus- 
ter S urvey, Postman and collaborators (Postman et al. 



1996 , P96 hereafter) wrote a matched filter algorithm to 



identify cluster candidates by using positional and pho- 
tometric data simultaneously. Slight variants have been 
proposed hence, with differences lying in some details of 
the detection processes (Kawasaki et al. 1998 ) or in the 
generalization of the algorithm so as to render it applica- 
ble to any type of survey data (Kepner et al. 1999| ). Our 
work, instead, tries to improve on the P96 algorithm in 
the sense of removing the a priori assumptions that are 
implicit in the Postman et al. technique. 

The P96 algorithm relies on the choice of both a given 
cluster profile - modificd-Hubble or King-like (hereafter 
often referred to as King filter for simplicity) - with a 
typical cluster scale - the core radius r c - and a cut-off 
radius, and of a typical M* - the chosen parameterization 
of the Schechter function. Both quantities, r c and M* , are 
rigidly coupled to detect cluster candidates and assign 
them a tentative redshift. Each one of these quantities 
being a function of distance (and thus cosmology), their 
implementation on any algorithm implies a dependence on 
the adopted cosmology and on the chosen amount of evo- 
lution, not to mention the dependence on the particular 
values actually chosen for r c and M*. 

In our new algorithm the spatial and luminosity 
part of the filter are run separately and succes- 
sively on the catalogue, with no assumption on the 
typical size or typical M* for clusters. In fact, these pa- 
rameters intervene in our algorithm only as a typical an- 
gular scale and typical apparent magnitude m*, bearing 
no ties to fixed physical scales nor to absolute magnitudes 
through redshift dependence. This has the consequence of 
removing the need for a choice of fixed physical values 
for these two quantities, and for a choice of evolutionary 
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models. Moreover, the fact of not coupling the space and 
the luminosity parameters also enables us to reach higher 
values of completeness: a candidate can be retrieved even 
if it does not flag a maximum likelihood at the very same 
redshift value simultaneously for both the space and the 
luminosity distributions (a situation that would lower its 
global likelihood when using the P96 algorithm). 

One further advantage consists in the local estimate 
of the background for each cluster candidate. In the 
approach of P96, the magnitude distribution is assumed 
constant all over the catalogue and so is the background 
spatial density. This may be hazardous especially in what 
concerns the spatial part, where local inhomogeneities 
may hamper the detection in shallower regions of the 
catalogue. Due to varying observing conditions, a part 
of the EIS data, at least, is reported to be non uniform 



(Olsen et al. 1999b), so this feature of our algorithm is 



of great help, as it allows us to tackle the problems of 
inhomogeneous data sets and to achieve higher levels of 
completeness. 

In our algorithm local enhancements in the projected 
galaxy density are first selected through the Gaussian 
filtering, and a "spatial probability" of them being spuri- 
ous is computed. Subsequently, the maximum-likelihood 
"filtering" on the apparent magnitudes searches for 
the presence of a Schechter distribution superimposed 
on the local background (which follows a power-law 
instead). This second step leads to the assessment of a 
corresponding "luminosity probability" that, multiplied 
by the spatial one, produces the final probability of each 
candidate being a spurious one. This means that the 
lower this value, the more confident the candidate. Notice 
also that this final quantity is always lower than either 
one of the partial spatial or luminosity probabilities as it 
is their product. 

We will now discuss in detail the components of the 
algorithm. 

2. 1 . The spatial filter 

In the spatial part we have chosen to work with a Gaus- 
sian filter. This choice, while avoiding to specify too 
particular a shape for the profile of the cluster candi- 
dates, brings along all the advantages and favorable math- 
ematical properties of the Gauss function. It is worth re- 
membering that the gaussian function in Fourier space is 
more compact than the King function. As a consequence 
the convolution with a random distribution of galaxies (a 
white noise in fourier space) will produce less spurious 
detections for the Gauss filter than for the King filter. 

Regarding completeness, it can be interesting to show, 
using simulations, the relative advantages of the Gaussian 
filter with respect to the King filter. A simple way of do- 
ing it is to simulate an area with randomly distributed 



field galaxies where to embed clusters, and apply to this 
area the two spatial filters in question to perform cluster 
detection and to compare directly completeness for both 
of them. 

Comparing the two filters is not straightforward, 
though: their shapes are intrinsically different so that 
there is no direct correspondence between the respective 
typical scales (r c and <7 a ng)- The King-like profile has a 
slimmer central peak with broad wings while the Gaus- 
sian shows the well known bell-shaped form. However, 
when cross-correlating the two filters, the highest signal is 
achieved if we consider a ang ~ 2.52 x r c , suggesting that 
a factor of 2.52 should be used in the simulations when 
comparing the relative efficiency of the Gauss versus King 
filters. Besides, since when searching for a spatial density 
enhancement one has to move in steps of dimension com- 
parable to the typical scale of the filter, this also implies 
that the Gaussian filtering requires a factor of (2.52) 2 ~ 6 
less search points with respect to the King filtering. 

As we have already mentioned, in our search we will 
not fix any typical physical size for our filter, but we will 
use, instead, a range of angular sizes u ang , namely from 
~ 0.35 up to ~ 1.42 arcmin in five steps of ratio 
These values were chosen bearing in mind the range of 
reasonable dimensions of the cluster candidates spanning 
the redshift interval that we expect to probe with the EIS 
data (z ^ 1.2). The corresponding five r c values to be 
used in the spatial search with the King-like filter are 
those obtained using the relationship a ang ~ 2.52r c . 



2.1.1. Simulations 

The simulations performed for this comparison were done 
as a function of the signal-to-noise threshold chosen, of 
the cluster profile, richness and redshift. 

As the final goal was to apply our algorithm to the 
EIS data, to simulate the galaxy background distribution 
we used the data from the EIS itself (see section ||), by 
random shuffling galaxy positions and magnitudes within 
the limits of the survey. 

The contamination rates were assessed by running the 
Gaussian filter and the King filter on these pure back- 
ground fields. 

In what regards completeness, this rate was assessed 
by running the two filters on the background fields plus 
cluster galaxies. This means that, for each simulation, a 
cluster of galaxies with different characteristics was placed 
randomly within each background frame, with the only 
constraint of being at adequate distance from the frame 
borders (ie. a distance larger than both 10 times the <J ang 
used for the search as well as 5 times the angular size of 
the embedded cluster), in order to avoid border effects. 
The characteristics of the clusters were defined according 
to the following prescription: 
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The cluster mean redshift: ranging from z — 0.2 till 
z = 1.2 at intervals of Az = 0.2. 

The surface density profile: a power law in radius of the 
form r 13 . The index j3 was set to three different values 
(0, —1 and —2) when testing the spatial filter only (as 
is the case of this section) . But it was fixed to the typ- 
ical in-between case (3 = —1.4 for the overall simula- 



tions (further ahead, in section 2.5, when we will apply 



the complete algorithm - that is, both spatial and lumi- 
nosity filters - to the mock catalogues); this provided 
a compromise that also allowed a direct comparison 
with the P96 results. Furthermore, we introduced a 
small central smoothing region in the profile, with the 
typical size of a cluster cD galaxy (r smoot h ~ 35 kpc, 
the average diameter of the two central giant galaxies 
in Coma, as provided by the NED-database) , to avoid 
a cusp. Both r smoot h and the total radial size of the 
cluster, 1 Mpc, were translated into the corresponding 
angular sizes according to the redshift chosen in item 
1. Doubling the cluster size to 2 Mpc did not change 
the results of our simulations in terms of completeness. 
The luminosity function - a Schechter, which is gen- 
erally adopted to describe th e lum inosity distribution 
of cluster galaxies (Schechter 1976| ), but with parame- 
ters determined by Colless ( 1989| ) for a set of 14 clus- 
ters of mean redshift < z >~ 0.0851 observed in the 
Bj band. Thus, the faint end slope is a = —1.25 and 
the characteristic apparent magnitude, at that mean 
redshift, is m* r — 14.8 when converted to the I-band 



(Fukugita et al. 1995). The m*j used in each simula- 
tion was changed according to the redshift (item 1) 
and to the k-correction applied (see item 5 below). 
The richness, as given by n* , the number of galaxies 
brighter than the characteristic magnitude m* . Three 
cases were considered: n* — 30, 50 and 80 (following 



Schechter 1976 ). These can be roughly ident ified w ith 



1958 



see 



Abell richne ss cla sses R— 1,2 and 3 (Abell 
also BahcaU |l98f). 
The dominant morphological type, as implied by the 
k-correction that is used. Actually only one extreme 
case was tested: clusters entirely composed of elliptical 
galaxies. We deliberately chose to ignore both the evo- 
lutionary corrections and the possibility of having later 
morphological types as the cluster dominant popula- 
tion, since this would correspond to a more favorable 
scenario of rendering cluster members brighter, which 
would only facilitate their detection, k-corrections have 
been provided by the spectro-photometric evolution- 
ary model PEGASE of Fioc & Rocca-Volmerange 



(1997) tuned to the cosmology parameters adopted. 



Once these characteristics are fixed, we can compute, 
for each cluster, the number of member galaxies that are 
observable within the apparent magnitude limit. Positions 
and magnitudes are next assigned to these cluster galaxies 
by random shots drawn from items 2 and 3, respectively. 



This was done in sets of 100 simulations and the two 
spatial filters were applied to these mock catalogues. Both 
field and cluster galaxies were limited to a magnitude 
cut-off of 1=22.0, the one we adopted also for the EIS 
data. 

In terms of completeness, the results are summarized 
in figures [l] and ||. 
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Fig. 1. Comparative performance, in completeness rate, 
between the spatial filter used in our new algorithm with 
the five adopted angular widths (u ang ) - filled squares - 
and the P96 King-like spatial filter with the corresponding 
five core radii - empty circles - such that r c = l/2.5a ang . 
Each point represents the mean value of a set of 100 boot- 
strap simulations. Clusters are elliptical dominated, have 
radial profile and richness class R ~ 1 (left), R ~ 2 
(center), and R ~ 3 (right). See text for further details. 

According to the plots in figures and |, the com- 
parative performance in completeness rate shows that the 
Gaussian filter is more efficient in most of the cases than 
the King filter, and this is especially true in particular in 
the difficult case when the cluster profile is less peaked. 

2.2. The luminosity filter 

To assess the signal given by the luminosity information 
we apply a maximum-likelihood technique on the appar- 
ent magnitude distribution of the candidates previously 
found with the spatial filter. We use as functional form 



for the luminosity function a Schechter ( 1976 ) function 
expressed in apparent magnitudes, thus avoiding the P96 
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Fig. 2. Same as figure | but now the richness class has 
been fixed to R — 2 and different r^ radial profiles are 
compared in terms of filter completeness performance. 
Left panels represent (3 = clusters, in the middle stand 
j3 = — 1 clusters and the right panels show results for 
(3 = —2, steeper profile, clusters. 



choice of an intrinsic M* and the choice of a given cosmol- 
ogy (H and q ) and morphological content (translated in 
the assumption of a given k-correction and, eventually, of 
an e-correction as well). Also, unlike P96, we have decided 
to implement the generalized treatment for th is likelihood 
filter, following Schuecker & Bohringer ( |l998[) . This com- 
plete treatment uses the exact mathematical equations, 
allowing to keep track of all possible errors that would 
affect any eventual direct redshift estimate (like the one 
performed by P96). It does not add major computational 
effort nor time. Besides, it also renders unnecessary the 
final "cluster signal correction" demanded by the complex 
approximations of the P96 procedure. 

W e thus took equation (8) of Schuecker & Bohringer 
( |l998|) and adapted it to compute the likelihood for this 
luminosity part of our algorithm. Notice that, as we de- 
couple it from the spatial part, the galaxy distribution 
profile is suppressed from the original formula so that the 
final log- likelihood ratio, for a given m* , is: 



hi 



^S+bkg 



L 



bkg 



N 



= E ln 



N - A 

N 



A (him. 



(N - A) b(mi) 
(1) 



where L stands for the likelihood function, S and 4'{ rri i — 
to*) refer to the Schechter function parameterized by a 
given m* , bkg stands for the background galaxies and b(m) 
for the differential magnitude number counts of the back- 
ground galaxies; b(m) is mathematically described by the 
fit to the magnitude distribution of all galaxies in the cat- 
alogue. N is the total number of cluster and field galaxies 
present inside a 2.5a ang radius circle, where we also esti- 
mate a coa rse richness parameter A (as will be detailed in 
section 2.3). 

Note that d> has been normalized, ie. : 



(to — to*) dm 



1 



(2) 



mu m being the limiting magnitude of the catalogue. 

For each candidate we compute the log-likelihood ratio 
(equation ^) for different to* values, set within the lim- 
its of magnitude of the galaxy catalogue by steps of 0.1 
magnitude. The slope of the Schechter function is fixed 
to a = —1.25, according to the mean value derived by 
Colless ( 1989 ) and in agreement with typical cluster lu- 
minosity function parameters found in the literature since 
the work of Schechter (1976). By identifying the absolute 
maximum in the distribution of the log-likelihood ratios 
relative to all to*, we automatically select the output to* 
for that candidate. 



2.3. Detailed description of the algorithm steps 

In this section we present a detailed description of the way 
the algorithm works when applied to the data. 

We note that we shall often interchange probability 
with the complementary Gaussian percentile and the as- 
sociated level of standard deviations, referred to as S/N 
throughout this paper. This approximation was used also 
in P96 so, by adopting it as well, we render comparisons 
among our results easier. 

The algorithm starts off by reading the coordinates 
and magnitudes of the galaxies present in a magnitude 
limited catalogue. Then, for each one of the five selected 
angular apertures a ang of the Gaussian, the code: 

1. Builds the search grid for the spatial filter analysis. 
The spacing between adjacent grid points is a ang but 
no grid point is allowed within 3cr ang of the catalogue's 
borders in right ascension nor in declination. 

2. Computes a first hand estimate of the local background 
density associated to each grid point by counting the 
number of galaxies in a circle with a 6a ang radius cen- 
tered at the grid point in question. (All cases of grid 
points distant less than 6a ang from the catalogue's bor- 
ders are obviously corrected by shifting the circle.) No- 
tice that, in this way, the background spatial density is 
not assumed constant all over the catalogue, but is es- 
timated around each point. This characteristic proves 
to be particularly efficient in dealing with data - like 
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those of EIS in at least some of the patches - where 
varying observing conditions result in non homogene- 
ity of the galaxy catalogue (Olsen et 



final global probability, P to t , for the detection of each 
cluster candidate. This is simply done by multiplying 



1999b) 



3. Performs the Gaussian filtering at each grid point con- 
sidering all galaxies within 3a ang from it and weigh- 
ing their relative position through the Gauss function. 
The resulting output signal is normalized (respective 
mean subtracted and division by the corresponding 
standard deviation), obtaining the S/N appropriate 
for each grid point. 

4. Selects all grid points which flag a local maximum, 
it. , for which the spatial signal is greater than that 
of all their immediate grid neighbours. No points at 
the extremes of the grid are allowed in this choice, 
thus avoiding "border effects". The signal must also 
be above a given threshold, imposed as an input in the 
programme. We set this spatial detection threshold to 
2.5 (the equivalent to the 98th percentile). 

5. For these maxima, our algorithm proceeds to the fine 
centrage: the Gaussian filtering is redone, maintaining 
the filter angular size, but building a finer grid in an it- 
erative procedure (down till an eighth of <J ang ) around 
the tentative candidate. The fine tuned center is fi- 
nally chosen as the position where the spatial signal is 
maximized in the ensemble of the points of this denser 
grid. 

6. The five different scale catalogues are cross-correlated 
in order to remove double detections. Detections are 
considered double if the distance between their centers 
is lower than or equal to the mean of their scales. In 
this case, only the highest signal detection is kept. 

7. For each selected candidate, the next step is to obtain 
a preliminary estimate of its richness, A, to be used in 
the computations of the likelihood of the luminosity 
part (see equation Q). A is the approximate number 
of galaxies statistically belonging to the cluster can- 
didate. We obtain it by subtracting to the N galaxies 
inside a 2.5a an g radius the local background, estimated 
in a annulus around each candidate. This ring has an 
inner radius of 5<7 ang and its area totals 0.06 deg 2 . 
These values were chosen after evaluating the galaxy 
density distribution around several candidates in the 
final EIS catalogue. 

8. For each candidate we compute the log-likelihood ratio 
(equation [l]) for the set of to* values and determine 
the one that maximizes it. The probability associated 
to this value is obtained by bootstrap: the catalogue 
is randomized, the spatial filter is run on it and the 
distribution of the log-likelihood for the flagged spatial 
detections is computed thus allowing an estimate of the 
probability of having an L max greater or equal than 
that measured for our candidate. 

9. At this point we possess the a ang and to* that inde- 
pendently maximize both the spatial and the luminos- 
ity signals and their respective significance, P space and 
Plum- We can combine these probabilities to obtain a 
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Fig. 3. Results from the simulations (that allowed us to 
assess the overall completeness rate - see section 2.5) show- 



ing the comparison between the redshift value of the sim- 
ulated clusters and the output to* obtained from our al- 
gorithm. The curve shows the relationship between z and 
to* used to define the input to* in our simulations (see 
text for further explanations). Each point is the mean of 
100 simulations, la error bars are drawn. 



2.4- Richness estimate 

The estimate of the richness of our cluster candidates is an 
important step, as we would like to quantify the percent- 
age of groups versus rich clusters that we will be detecting. 

Whichever richness estimate one decides to adopt, a 
necessary step is the determination of a physical radius 
for each cluster candidate which, in turn, requires the es- 
timate of its redshift. Unlike P96, our method does not 
provide a direct estimate of the redshift for our candi- 
dates. We obtained such estimate using to* and its rela- 
tionship to M* , whose adopted value is that reported by 
Colless ( |1989D , to be redshifted and corrected using the 
k-corrections typical of ellipticals, as determined by Fioc 
& Rocca-Volmerange (1997). Simulations show that the 
value of to* obtained through our algorithm retrieves well 
the M* introduced in the simulations - see figure ||. On 
the contrary, the value of a ang retrieved by our algorithm 
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does not correlate well with redshift: the angular size cho- 
sen by the algorithm has a large scatter with respect to 
the physical size of the clusters entered in the simulations. 

Having thus obtained a rough estimate of the candi- 
date's redshift, we can also estimate its rich ness, adopting 
the iV .5 parameter proposed by Bahcall fll.981 ). It con- 
sists in counting the number of cluster member galaxies 
brighter than 7713 + 2 (where is the magnitude of the 
third brightest cluster galaxy), located within a projected 
radius of 0.5 Mpc from the cluster center, the typical size 
of a cluster X-ray emitting region. Background correc- 
tion is estimated from similar counts performed in the 
entire catalogue region, using 0.2 magnitude bins. iVo.s 
was found to be well correlated with the system's veloc- 
ity dispersion v r by Bahcall (1981, see also eg. Alonso et 
al. |1999 for further studies on this type of correlation) 
for a sample of 26 nearby clusters ranging from rich sys- 
tems (v r ~ 1500fcms _1 ) to small groups of galaxies (v r ~ 
lOOfcms" 1 ) and has already been used in the literature for 
richness estimates of distant clusters (Deltorn et al. 1997). 
It has the advantage of avoiding the higher uncertainties in 
background subtraction present when adopting the stan- 
dard Abell ( |1958[ ) estimate, which is performed within a 
radius of 3 Mpc (see eg. P96 for comments on this). In 
fact, is relatively insensitive to uncertainties in ei- 

ther the background correc tion o r in the exact position of 
the cluster center (Bahcall 1981). iVo.5 can be translated 
into the typical Abell rich ness c lasses using the indicative 
relation given by Bahcall ( 198l| ), NAbeii — 3.3 x iVo.g, and 
then checking table 1 of Bahcall ( |1988[ ) . This means that 
eg. No. 5 < 15 corresponds to Abell richness class R ~ 0, 
while N0.5 > 24 is equivalent to R <; 2. 



2.5. Expected final overall rates of completeness and 
contamination - simulations 

We ran a new set of extensive simulations in order to 
assess the overall efficiency of the algorithm. The layout 
of these simulations being the same as the ones that were 
used for estimating the completeness and cont amin ation 
rates of the spatial filter alone (see section [2.1.l| ), we 
shall not repeat their details. Here, it will be enough to 
mention that the mock catalogues (pure field for contami- 
nation rates or field plus cluster galaxies for completeness 
purposes) were now submitted to the algorithm, a nd not 
only to its spatial filter as had been done in section 2.1.1 



The results regarding the overall completeness rates 
are shown in figure |^. Even setting the detection threshold 
at S/N ~ 4.0 - which renders contamination rates negli- 
gible (see below for the final overall rate of contamination 
directly assessed on the EIS data) - we are able to achieve 
a completeness of ~ 95% until z ~ 0.9 for richness 2 
Coma-like clusters, and of ~ 97% up to z ~ 0.8 for 
poorer richness class 1 systems. For R = 3 clusters such 
high values of completeness hold even in the redshift bin 



z ~ 1.0, always for the no-evolution case. These results 
improve on the values obtained using the P96 algorithm 
(see their figure 20). In fact, while we set the Schechter 
slope of the simulated clusters to a steeper value (—1.25 
instead of —1.1), the fact of having adopted a brighter 
magnitude limit (22.0 versus P96's 22.5) renders a direct 
comparison of both works plausible. We checked this by 
integrating both Schechters for the two magnitude limits, 
thus getting 38(us) versus 50(P96) galaxies for R=l 
clusters, 63 versus 82 galaxies for R=2 clusters, and 100 
versus 131 galaxies for R=3 clusters. Moreover, do notice 
that the cut-off S/N is actually very different: while we 
set this threshold to 4cr, P96 report as detection limit the 
95th percentile [ie. approximately the 2a level) in their 
simulations. 
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Fig. 4. Percentage of completeness for Abell richness 
classes R = 1 (poor), R = 2 (intermediate) and R = 3 
(rich) clusters as established through simulations. Clus- 
ters have a power law radial profile of type r 
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arc 



dominated by ellipticals, and have a typical luminosity 
function as determined by Colless (1989). We applied el- 
liptical k-corrections (Fioc & Rocca-Volmerange 1997) to 
cluster galaxies. 



As for the expected contamination rate, we ran our 
algorithm on the EIS catalogue after having randomized 
galaxy positions and magnitudes. This time, no clusters 
were added to the galaxy background distribution, as we 
were interested in the false detection rate associated to a 
given threshold S/N. These simulations supply ~ 14 spu- 
rious candidates per square degree if the threshold is set 
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at S/N — 3, and ~ 1.3 per square degree if it is set at 
S/N = 4. P96 report an estimated contamination rate for 
their final catalogue (of 79 candidates) of at most 30% 
in their 5.1 square degree area, corresponding to a sur- 
face density of contaminants of ~ 5 per square degree for 
PDCS. However, a direct comparison between the two re- 
sults is very difficult since the data are different and the 
algorithms have very different selection methods. 



Finally, and for technical reasons, we trimmed each 
patch so as to have clear-cut rectangular edges. This trans- 
lates into the limits given in table [l], and in the follow- 
ing number of galaxies: 57366 for patch A-even, 57553 
for patch A-odd, 25456 for patch B-even and 25445 for 
patch B-odd. In summary, we have four catalogues, two 
(the even and the odd) for each patch, where we shall 
independently run our algorithm. 



3. The catalogues 

Our final catalogues contain the following entries for each 
one of the candidate clusters found in ElS-wide patches 
A and B: right ascension and declination (equinox 2000), 
indicative angular size 2.5a ang (in arcmin), combined 
total probability expressed as a S/N, to*, redshift and 
richness (iVo.5) estimates. 
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Table 1. Limits of the regions to be analyzed by our al- 
gorithm (coordinates are expressed in equinox J2000). 



3.1. Defining our data set 

ESO provides on the web the single frame catalogues 
for the even and the odd coverage of the EIS patches 
( http: / / www.eso.org/ science /eis / ) . We thus had to assem- 
ble these data and build two galaxy catalogues (henceforth 
the even and the odd) for each entire patch, that were the 
basis for the cluster searches. 

After careful inspection of the EIS data character- 
istics, we decided to adopt the following parameters to 
define the galaxy catalogue on which to perform our 
automated search. The star/galaxy separation line was 
drawn at 0.9: only the objects with stcllarity index lower 
than this value were kept, irrespective of their magnitude. 
Let us notice that this criterion slightly differs from 
that adopted by Olsen et al. (1999a, 1999b), where all 
the objects at magnitude fainter than 21 were kept, 
irrespective of their stellar index, introducing a sharp 
step at 1=21 in the galaxy number counts. 



We excluded objects with / > 4 and nflag/npix > 0.01 
(where / are the SExtractor flags, npix the number of 
pixels above the analysis threshold and nflag the number 
of pixels flagged by WeightWatcher - see Nonino et al. 



1999 for details). The layout of the EIS mosaic is such 



that adjacent odd/even frames have a small overlap at the 
edges, especially in declination and we further did cure 
double detections in frame borders (ie. for pixel values 
X < 100, X > 1950, Y < 300 or Y > 1600) using a 2 
arcsec separation cut-off. 

We limited the data set to 18 < I < 22.0, aiming at 
guaranteeing a high level of completeness of the data. 



We note that Nonino et al. (1999) estimate that the 



single-frame odd and even catalogues are 80% complete 
to I = 23 for a typical EIS frame. 



3.2. Defining our S/N thresholds 

The choice of a threshold value in S/N for cluster detec- 
tion has to be a compromise between the will to achieve 
a high level of completeness and the need to avoid unac- 
ceptable high values of contamination. 



In section 2.5 we have shown that a cut-off value of 
S/N = 4tr gives a number of contaminants of ~ 1.3 per 
square degree, while guaranteeing a high level of complete- 
ness (> 97% for all three richness classes R = 1 — 3 up to 
a redshift of 0.8 - see figure ^). 

If we had in our hands two good quality even and odd 
catalogues for each patch, we could just apply such cut- 
off on both catalogues and get two final lists of cluster 
candidates that should differ by a small number, caused by 
false detections in regions where eg. SExtractor detected 
spurious objects in the spikes of bright saturated stars or 
satellite tracks. 

Unluckily, this does not appear to be the case for the 
data, especially in what concerns patch A, whose CCD 



frame quality is quite variable (Olsen et al. |1999q ). As a 
result, applying a blind cut-off of 4a to both catalogues 
would be too restrictive a choice, and therefore we decided 
to follow a more flexible approach. 

For each patch we first built, for the corresponding 
even and odd galaxy catalogue, a catalogue of cluster can- 
didates with a S/N cut-off of 3. Then, from these cata- 
logues we selected: 

Class 1 candidates: present in both even and odd 
catalogues (which implicitly means that they have 
S/N > 3.0cr in both catalogues) with S/N > 4. Oct in 
at least one of them; 

Class 2 candidates: present in one catalogue only but 
having S/N > 4.0ct. 







In this way, and for each patch, we do not reject those 
candidates that are present only in one catalogue; in fact, 
very often, these single detections include, other than fla- 
grant false candidates, systems non detected in one of the 
catalogues due to shallower images, worse seeing condi- 
tions or general lower photometric quality. 

In terms of contamination, our choice is equivalent 
to applying a S/N > 4.0 cut-off to the even/odd cata- 
logues separately and the same holds for completeness, but 
with the further constraint for class 1 candidates of being 
present also in the odd/even catalogue with S/N > 3.0 
(no further constraint for class 2 candidates). 

3.3. Cluster candidates 
3.3.1. Patch A 

Class 1 candidates in patch A total 41 (see table ^), which 
translates into a projected number density S ~ 13 clus- 
ters per square degree (this patch covers ~ 3 square de- 
grees). Three further class 1 candidates were detected in 
patch A but are not listed in Table ||. Two were cases of 
a bright star split into many objects and causing false de- 
tections, and a third one was a nearby low surface bright- 
ness galaxy again split into many components. The total 
expecte d n umber of spurious detections in this area (see 
section 2.5), is < 4. 



Table |2] lists, for each candidate, its right ascension, 
declination (both at the year 2000), a typical size (2.5a an g, 
where cr ang is the size of the spatial detection filter in ar- 
cmin), the S/N in the even and odd galaxy catalogues, the 
output to* , the estimated clust er rcdshift and the richness 
indicator Nq.s (Bahcall, 1981). The last column reports 



identifications with other catalogues, whenever is the case. 

The indicated size is the mean value of what was mea- 
sured from the even and odd catalogues respectively. The 
same was computed for all the other quantities in the ta- 
ble (S/N excluded, obviously), except whenever at least 
one of the two next cases occurred: (1) The redshift es- 
timates from the even and odd catalogues differed in ab- 
solute value by more than 0.1. This happened for a small 
percentage of the entries (8 over 41 ie. 20 %), suggesting 
that our 'a posteriori ' redshift estimate is quite robust "in- 
ternally" (between even and odd catalogues). (2) The rich- 
ness estimate was hampered because the candidate was 
attributed to 3 + 2 > 22 (the adopted magnitude limit) in 
the even catalogue (noted (e)) or in the odd (marked (o)) 
or in both ((eo)). In these cases, results issued from both 
even and odd catalogues were listed. 

Furthermore, candidates marked with an asterisk have 
only a lower limit for the S/N, obtained using the spatial 
information alone. For these candidates we preferred to 
avoid quoting the total S/N as their Pi um was extremely 
low, possibly an artifact due to bad data quality: a 
careful inspection revealed that local faint galaxy counts 
were significantly below those of patch A as a whole 



thus causing, as a consequence, an overestimate of the 
significance of the bright galaxy excess. 

As for class 2 candidates in this patch, table [| lists 
29 detections: those that survived a visual scrutiny aimed 
at rejecting obvious cases of false signals (satellite trails, 
bright star refraction spikes and other artifacts). For each 
cluster candidate the same quantities as in Table || are 
listed. 

The reason for a candidate - not being an obvious ar- 
tifact of the data - to be present only in one of the two 
catalogues of patch A has to be related to the worse quality 
of the corresponding CCD frames in the catalogue where 
it is missing. The variables at work are numerous (such 
as seeing conditions, background noise, etc.) and not al- 
ways easy to define. Often, a noisy background results in 
a higher number of - possibly spurious - galaxy detections 
by SExtractor, while a worse seeing invariably results in 
a lower number of galaxy detections. 

The large number of class 2 detections is an indication 
of the large range in quality of the CCD data available for 
patch A. In fact, if we check the ratio of class 2 to class 1 
objects, we get 29/41 or, in other words, class 2 candidates 
represent a percentage of ~ 41% of the total sample. Olsen 
et al. report, for the same ratio, 15/20, which translates 
into ~ 43% of their whole sample. 

Figure |5| shows the sky distribution of our class 1 and 
class 2 candidates, circles and squares respectively. 

Adding class 1 to class 2 candidates provides a total 
of 70 cluster candidates for patch A, corresponding to a 
surface density of 23 per square degree, with 1.3 expected 
spurious detections per square degree. Out of this total, 
40% can be identified as richness class objects and the 
same percentage encloses R ~ 1 systems, while the re- 
maining 20% seem to be richness R ~ 2 clusters - see the 
upper panel of figure |l The corresponding redshift dis- 
tribution for the total Patch A sample is plotted in the 
upper part of figure [?]: the distribution of class 1 candi- 
dates does not seem to differ significantly from that of 
class 2 candidates, both covering the redshift range 0.3 to 
1.0 (the excess tail present for class 2 candidates in the 
last redshift bin may be an artifact, due to those candi- 
dates whose to* approaches (or would exceed) I = 22, the 
limiting magnitude of our galaxy catalogue). 

A detailed comparison of our catalogue of cluster can- 
didates with that of Olsen et al. (1999b) will be presented 
in section ^. Here, let us notice that in patch A there are 
four cluster candidates detected by previous independent 
works and listed in NED. Abell cluster ACO1055, that 
was selected as a class 1 candidate by us with a high S/N. 
C12245-4002 (LP96) of Lidman & Peterson ( |i996| ) is well 
identified with our entry 22 ,l 47 m 54.2 s , -39°46'33.7". Two 
EDCC cluster candi dates (EDCC163 and EDCC169; 
Lumsden et al. 1992Q do not turn up in our sample. We 
decided to point our algorithm directly at the NED co- 
ordinates of these two clusters while using the galaxy 
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Fig. 5. Map of patch A candidates. Circles indicate class 1 systems while squares designate class 2 candidates. 
Coordinates and sizes (proportional to a ang ) are mean values from even and odd catalogues. Do notice that some 
class ,8 candidates (namely at, roughly, (339.3°, —40°) and (340°, —40.35°)) may misleadingly appear to have escaped 
our double detection elimination procedure but the reason is that apparently superimposed detections are actually 
unmatched (class 2) ones flagged in different catalogues (even or odd). 



catalogue at Iu m — 23, just to check if their absence 
from our list could be due to the magnitude cut-off cho- 
sen. This approach managed to retrieve EDCCW9 with 
S/N even = 3.29 and S/N odd = 3.84, while EDCC163 is 
still missing. According to the Lumsden et al. ( 1992|) table, 
both candidates are considered to be poor R < systems, 
which may account for our results. 

Out of these 4 objects, only Abell ACO1055 has spec- 
troscopic confirmation (at z = 0.0322) and, consequently, 
a reliable richness estimate (R = 0) in the literature (cf 
NED database). The redshift we estimate for this system 
is higher (z even = 0.6, z a dd = 0.5) and the richness we 
compute (R ~ 2) is, consequently, considerably incorrect. 
Notice though that visual inspection indicates that a lot 
of spurious "galaxies" were detected by SExtractor in the 
spiral arms of the foreground bright galaxy ESO 345-G046 
(the effect is more serious for the even catalogue) , and this 
has possibly affected our results, namely the determina- 
tion of TO*. 



3.3.2. Patch B 

Class 1 candidates in patch B total 21, which translates 
into £ ~ 15 clusters per square degree (this patch cov- 
ers, approximately, 1.4 square degrees). This is a slightly 
higher density than that found for patch A, possibly be- 
cause of better data quality (patch B had better observ- 
ing conditions). Note that Olsen and collaborators also 
find some difference: they report 14 clusters per square 
degree in patch A while the value rises to 17.2 for patch 
B. This difference between both patches is, however, neg- 
ligible within uncertainties. 



Table [| lists, with the same format already used for 
patch A candidates, our class 1 patch B cluster candidates. 



The simulations of section 2.5 allow us to predict 
1.3 x 1.4 ~ 2 spurious detections associated to the 21 
candidates. 

Very much like for patch A, in 71 % of the cases (15 
out of 21), the redshift estimates from the even and odd 
catalogues differ in absolute value of less than 0.1. 

As for class 2 candidates in this patch, table |s| lists all 
8 detections that survived accurate visual inspection. 

Figure |^ shows the sky distribution of our class 1 and 
class 2 candidates, respectively. 

Also for patch B holds the discussion already made for 
class 2 candidates in patch A. It is to be noticed that now 
the ratio of class 2 to class 1 candidates is much lower than 
the one obtained for patch A: 8/21 or a percentage of 28% 
of the total sample (versus the previous 41%), suggesting 
that the quality of the CCD data available for patch B is 
indeed much more homogeneous. We note here that the 
same type of ratio produced with Olsen et al. 's candidates 
for this patch is 9/10 or, in terms of percentage of their 
whole sample, ~ 47%, which remains quite similar to the 
figure they produce for patch A. 

Adding class 1 to class 2 candidates gives a total 
of 29 cluster candidates in patch B, corresponding to a 
surface density of ~ 21 per square degree, with, again, 
the already mentioned 1.3 expected spurious detections 
per square degree. Out of these 29 systems, 62 % have 
estimated richness class R ~ 0, and 38 % have R ~ 1, 
as illustrated in the lower panel of figure ||. Notice that, 
the area probed here being smaller than the one of Patch 
A, at the same limiting magnitude this translates into 
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Fig. 6. Same as figure ^ but for patch B candidates. 



a smaller sampled volume. So it is expected that the 
median typical richness of the detected systems should 
move to lower values. The corresponding total redshift 
distribution is plotted in the bottom of figure fj]. The 
median of the redshift distribution for class 1 candidates, 
patch A and B considered together, is z — 0.65, similar 
to that of Olsen et al. despite having chosen a brighter 
cut-off in magnitude for the galaxy catalogue (22 versus 
their 23), while it is slightly deeper than that of the P96 
sample PDCS (z = 0.4). 

There is only one cluster candidate in patch B detected 
by previous independent works and listed in NED. It is 
cluster AC084, that was selected as a class 1 candidate 
by us. We estimate 0.40 for the redshift of this candidate 
(and a corresponding richness class R ~ 1), while spec- 
troscopic measurements in the literature hold a value of 
z = 0.11 (cf the NED database) and R = 0. Also for patch 
B, a more detailed comparison between our catalogue of 



cluster candidates and that of Olsen et al. (1999a) will be 
presented in section |[ 

4. Comparison with Olsen et al. 's results 

We can now do a thorough comparison with the results 
that Olsen et al. ( |l999a| , [19991: ) obtained by applying the 
P96 algorithm to patches A and B of EIS. In this way we 
will also be able to assess directly the relative efficiency 
of the two algorithms (P96 and ours) on the same set of 
data. 

We will have to do a double check: first which of our 
candidates are new with respect to the Olsen et al. lists, 
and then which of the Olsen et al. candidates are missing 
from our lists. 

For the comparison to be fair, we have to remember 
three important points: Olsen and co-workers did their 
search on the total galaxy catalogue, ie. up to magnitude 
/ = 23, while we preferred to limit it to / = 22. They 
applied a S/N cut-off of 3, while we adopted a S/N cut-off 
of 4 in our final cluster list and, finally, they also imposed 
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Fig. 7. Redshift distribution for the cluster candidates 
flagged in patch A - upper panel - and patch B - 
lower panel. Hatched histograms specify matched (class!) 
detections among the global distribution (empty his- 
tograms) . In this figure the z plotted for each cluster can- 
didate is the mean value obtained from even and odd cat- 
alogues. 



a cut-off in richness to their candidates, while we apply 
no such sort of selection. 

Bearing this in mind, we proceed with the compar- 
isons. 

Regarding patch A, in our class 1 list we have 13 
candidates in common with Olsen et al. ( |1999b[ ), while 
in the class 2 list we have only 3 candidates in common. 
Out of our remaining candidates, 6 from class 1 and 9 
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Fig. 8. Distribution of the richness estimate A0.5 for the 
cluster candidates flagged in patch A - upper panel - 
and patch B - lower panel. Hatched histograms specify 
matched detections among the global distribution (empty 
histograms). The dotted lines mark the approximate lim- 
its between Abell richness classes and 1 (at Nq.s ~ 15) 
and Abell richness classes 1 and 2 (A0.5 ~ 24). The iVo.5 
plotted for each cluster candidate is the mean value ob- 
tained from even and odd catalogues. 



from class 2 are in an area of patch A that they decided 
to cut-out from their search (the upper right corner of 
their figure 2) because of "obvious incompleteness" (sic). 
We did search also in this area, since our algorithm, 
using local background values, is better suited for dealing 
with inhomogeneous data. We therefore are left with 
22 class 1 candidates and 17 class 2 candidates present 
in our list and missing in theirs. Some of these missing 
candidates might have been discarded by Olsen et al. 
(1999b) because of the richness cut-off they imposed, but 
it should be noticed that amongst them there are also 
relatively rich systems: 5 with estimated richness R ~ 2 
and 25 with richness R ~ 1 (conservatively counting only 
our iVo.5 secure values), whose absence is thus not easy 
to justify. 



In what concerns patch B, we only have 7 candidates 



in common with Olsen et al. (1999a), all in our class 1 
list. We are therefore left with 14 class 1 candidates and 
all 8 class 2 candidates present in our catalogues and 
missing in theirs. 



On the other hand, there are candidates present in the 
Olsen et al. lists that we failed to select, and we wanted to 
check if our missing them could be related to the different 
magnitude cut-off of the galaxy catalogue used for our 
search or to the different S/N cut-off adopted. 

To investigate this we decided to check which would 
have been the S/N obtained for these missing candidates 
by our algorithm on the galaxy catalogue limited to their 
deeper magnitude cut-off (I = 23.0). To do this we simply 
targeted directly the search to the positions of each one of 
their candidates. 

By doing so, we were able to retrieve all their patch 
A candidates with a S/N > 3 in at least the even or the 
odd catalogue (the same selection criterion they adopted) , 
with the only exception of four of them. One of those 
was commented by Olsen et al. as being a "doubtful case 
based on the visual inspection of the coadded image" and 
to another one our algorithm could only attribute a lower 
limit S/N (from the spatial part). 

In patch B, with the same procedure, we assigned all 
their candidates a S/N > 3a except for three (one of 
them being also a case of lower limit S/N obtained from 
the spatial filtering). 

The plot of our S/N against theirs, for all the can- 
didates we retrieved (common candidates), is shown in 
figure ||, both for patch A (filled symbols) and patch B 
(empty symbols) objects. 
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Fig. 9. S/N obtained with our algorithm for the Olsen 
et al. (1999a, 1999b) candidates versus their own S/N . 
Filled symbols stand for patch A objects while patch B 
candidates are noted by empty symbols. The S/Noi sen = 
S/N us line is also marked. 
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One can notice that, for the majority of the cases, 
our S/N is greater than that obtained by Olsen et al. 
as would be expected: by decoupling the spatial and 
the luminosity parts of the algorithm, a candidate can 
be retrieved even when it does not flag a maximum 
likelihood at the very same redshift value simultaneously 
for both distributions, a situation that would lower its 
global likelihood when using the P96 algorithm. 

Another interesting plot is the one showing, for the 
common candidates, our m* estimate versus the redshift 



estimate by Olsen et al. (1999a, 1999b). Figure 10 plots 



the mean value (from the even and the odd catalogues) 
of m* produced by our algorithm against the redshift es- 
timated by Olsen et al. for their candidates in patch A 
(filled circles) and in patch B (empty circle s). Th e line 
indicates the redshifted value of the Colless ( 1989 ) local 
mean m* , affected by k-correction typical of ell ipticals (as 
determined by Fioc & Rocca-Volmerange 1997 ). The gen- 
eral distribution of the data points does seem to follow the 
curve's trend. The larger redshift systematically predicted 
by the Colless relationship we used is just a byproduct of 
the brighter M* estimated by Colless (I* ~ -22.2 for 
Ho = 100 km s _1 Mpc~ 1 ) versus the one used in P96 and 
by Olsen et al. (-21.28, also for H = 100 km s^Mpc" 1 ). 
This fact may also be contributing to the discrepancies 
between our estimates and the literature measures in the 
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Fig. 10. Checking the possibility of correlation between 
our m* (average value from even and odd catalogues) and 
Olsen et al. 's redshift estimate for their candidates in 
patch A (filled symbols) and in patch B (empty symbols). 
The line is the relation deduced from Colless (1989). 



case of the Abell clusters mentioned in section 3.3 



5. Summary and Conclusions 

In this paper we presented a new algorithm for cluster 
detection that improves on the P96 one, mainly by 
avoiding the need to assume a typical physical size or 
a typical M* for galaxy clusters. These two parameters 
intervene in our algorithm only as typical angular scales 
and a typical apparent magnitude m* , bearing no ties to 
fixed physical scales nor to absolute magnitudes through 
redshift dependence (and chosen cosmological model). 
One further advantage of our algorithm with respect to 
P96's consists in the local estimate of the background 
for each cluster candidate, particularly useful in the 
spatial part where local inhomogeneities (due to varying 
conditions during data acquisition) may hamper cluster 
detection in shallower regions of the catalogue. 

We applied this new algorithm to the ElS-wide 
database and did a thorough comparison with the results 
obtained using P96 on the same kind of data. As the EIS 
observing conditions varied considerably throughout the 
granted nights, there was a considerable spread in the 
data-quality of different EIS frames (Nonino et al. 1999Q . 
Both the particular features of the algorithm presented 
in this paper, that somehow compensate for the lack of 
homogeneity of the data, and the conservative limiting 
magnitude adopted for cluster search try to minimize the 



problems inherent to the EIS data, allowing to achieve a 
higher completeness level than that obtained using P96 
on the same data. 

In fact, we notice that the distribution of our patch 
A candidates seems more homog eneous than the one 
presented by Olsen et al. ( 1999b ): compare our figure 
[| with their figure 4, where they have a region around 
a° ~ 341.5, 8° ~ —40 devoid of detections, not to mention 
the fact that they intentionally left out a part of the 
patch when running their algorithm to "overcome" data 
problems related to the lack of homogeneity. 

Regarding the cluster candidate surface density, a 
precise comparison with the values presented in the 
literature is not trivial, due to different magnitude limits 
adopted for the search and to the fact that we did 
not apply any a priori selection on cluster richness (as 
P96 for example). Just as a qualitative comparison, let 
us remember that the PDCS group (P96) reports the 
detection of 79 cluster candidates over an area of 5.1 
square degrees till I — 22.5, or a number density slightly 
above 15. This value is consistent with the observed 
X-ray logN — logS around fluxes > 6~ 15 erg cm~ 2 s -1 



(see figure 3 of Rosati 1988). In this paper, adding class 



1 to class 2 candidates both for patch A and B we obtain 
a surface density of ~ 22 clusters per square degree, 
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out of which ~ 54% have R = 1 or above, which agrees 
well, within the uncertainties, with the number quoted 
by the PDCS group, possibly suggesting a slightly larger 
surface density for our candidates, especially when taking 
into account the deeper magnitude limit adopted for the 
search by the PDCS group. 

We are doing multicolor (B, V, R and I) observations 
of the sample of cluster candidates presented in this 



Escalera E., MacGillivray H.T., 1995, A&A, 298, 1 
Escalera E., MacGillivray H.T., 1996, A&AS, 117, 519 
Fadda D., Slezak E., Bijaoui A., 1998, A&AS, 127, 335 
Fioc M., Rocca-Volmerange B., 1997, A&A, 326, 950 
Fukugita M. Shimasaku K., Ichikawa T., 1995, PASP, 107, 945 
Gal R.R., de Carvalho RR., Odewahn S.C, Djorgovski S.G., 

Margoniner V.E., 2000, AJ 119, 12 
Gio ia I.M., 1997, pr oceedings of "The Young Universe" or 



astro- ph/9712003 



paper lismg the 3.6-m and 1NTT ESU telescopes. Kedslutt 
confirmation ot three ot our high rcdshiit [z ~ U.boJ 
candidates has already been obtained using FORS at 
the VLT, confirming the efficiency of the strategy adopted. 
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27/22 
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-40 00 44.1 


2.13 


5.00 


5.81 
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22/25(eo) 


EIS 2241-4001 


22 41 33.7 
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1.25 


3.67 


4.30 
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0.40 
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1.33 
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Table 2. Class 1 cluster candidates of patch A, ordered in right ascension. Coordinates (a,(5)2ooo arL d approximate 
angular size 2.5<7 onff (expressed in arcmin) are the mean values of the quantities issued from the even and the odd 
catalogues. We report values of to*, z and iVo.5 for both the even and the odd catalogue if the discrepancy amongst 
at least one of the first two quantities was large or if to 3 + 2 > 22 (22 being our adopted magnitude limit) in the 
even catalogue, signalled by (e), or in the odd one, (o), or in both, (eo), making the estimate of iVo.5 quite uncertain. 
Some further remarks on individual entries: we note that a bright galaxy at the center in the line of sight of candidate 
22 h 38 m 49.3 s , — 39°34'16.8" may have caused some (4 at the most) spurious detections, but it is hard to say as the 
background noise is enhanced in that area. Candidate 22 43 17.4, —39 51 54.3 is close to a large foreground spiral 
(ESO 345-G046) which arms were broken into pieces by the SExtractor procedure, thus producing spurious detections 
(see text). Finally, candidate 22 ft 44 m 40.4 s , — 39°45'29.4" is near a bright star, but galaxy detections arc not affected 
by it (maybe 4 of the objects detected in the cluster periphery are doubtful). 
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Table 3. Class 2 patch A cluster candidates. Columns and notes follow the scheme of the ones in table 
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21 


40/22 


00 


1.0/1.1 


12/13(eo) 


EIS 0050-2941 


00 


50 


09.4 


-29 25 


12.4 


3 


02 


3.43 


4.29 


20 


30/22 


00 


0.7/1.1 


9/17(eo) 




00 


51 


47.6 


-29 29 


20.2 


3 


54 


5.19 


4.97 




21.05 




0.90 


8 




00 


51 


49.6 


-29 45 


20.1 





88 


4.21 


4.29 


20 


10/20 


50 


0.7/0.8 


ll/15(o) 




00 


51 


50.5 


-29 38 


05.3 





88 


4.76 


3.36 


21 


80/19 


90 


1.1/0.7 


13/4(eo) 




00 


52 


49.3 


-29 27 


49.5 


1 


07 


4.83 


4.18 


21 


80/21 


90 


1.1/1.1 


17/15(eo) 


EIS 0052-2927 


00 


53 


07.1 


-29 47 


40.7 


2 


65 


4.04 


3.43 


21 


50/20 


60 


1.0/0.8 


11/11(0) 




00 


53 


38.1 


-29 28 


41.8 





88 


3.00 


4.09 


21 


10/21 


50 


0.9/1.0 


12/9(eo) 





Table 4. Class 1 cluster candidates of patch B, ordered in right ascension. Quantities listed and notes used are 
equivalent to those reported in table |. Candidate Q0 h 51 m 50.5 8 , -29°38'05.3" is near a bright star (and 4 of the 
objects detected in the cluster periphery are doubtful). 



18 



Lobo el al. , ElSily looking for distant clusters of galaxies 



0:2000 


52000 


2.5 (Jang 


$ 1 -^V even 


S/N odd 


m* 


Zest 


iVo.5 


Comments 


00 44 45.0 


-29 38 17.1 


0.71 


4.18 




21.30 


1.0 


20(e) 




00 47 21.9 


-29 43 16.9 


1.41 




4.90 


19.70 


0.6 


11 




00 47 56.7 


-29 26 00.1 


0.35 




4.68 


22.00 


1.1 


9(o) 




00 48 13.5 


-29 44 44.3 


1.41 


> 4.21 * 




18.00 


0.3 


12 




00 48 18.2 


-29 49 09.5 


0.35 




4.27 


20.10 


0.7 


11 




00 48 22.1 


-29 45 21.4 


0.35 


5.25 




18.00 


0.3 


18 




00 52 00.5 


-29 50 58.3 


0.35 




4.00 


21.80 


1.1 


11(0) 




00 52 43.6 


-29 45 45.3 


0.35 


4.18 




21.70 


1.1 


14(e) 





Table 5. Class 2 patch B cluster candidates. Columns and notes are the same that were used in table 



